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Abstract. One aspect of supernova remnant evolution that is relatively unstudied is the influence of an AGN 
environment. A high density ambient medium and a nearby powerful continuum source will assist the cooling 
of shocked ejecta and swept-up gas. Motion of the surrounding medium relative to the remnant will also affect 
the remnant morphology. In an extension to previous work we have performed 2D hydrodynamical calculations 
of SNR evolution in an AGN environment, and have determined the evolutionary behaviour of cold gas in the 
remnant. The cold gas will contribute to the observed broad line emission in AGNs, and we present preliminary 
theoretical line profiles from our calculations. A more detailed comparison with observations will be performed in 
future work. The SNR-AGN interaction may be useful as a diagnostic of the AGN wind. 
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1. Introduction 

Active galaxies produce prodigious luminosities in tiny 
volumes, display very strong cosmological evolution, and 
are known for the diversity of their behaviour (e.g., 
Osterbrock & Matthews 1986). Their activity arises from 
the release of gravitational energy from accretion onto a 
supermassive (~ 10^ — 10^ M0) black hole and/or the cu- 
mulative effects of short-lived episodes of nuclear star for- 
mation, and is often accompanied by nuclear winds and 
jets. There is evidence (e.g., Williams & Perrv [TUM|l that 
intense star formation and nuclear activity are related, 
and starbursts are one possible way to fuel the central 
black hole. In turn, the nuclear wind and the extreme ra- 
diation field act back on the starburst components, in- 
fluencing stellar winds (and thus stellar evolution), wind 
blown bubbles, and supernova remnants. For the most lu- 
minous active galactic nuclei (AGNs) to be powered by ac- 
cretion requires that the density of the interstellar medium 
be enhanced, for example by galaxy encounters or mass 
loss from a dense nuclear stellar cluster. One of the most 
important questions concerning the AGN phenomenon is 
therefore the connection between starburst and nuclear 
activity. 

A defining characteristic of AGNs are their possession 
of strong (and often very broad) line emission. This pro- 
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vides detailed information on the physical conditions right 
down into the AGN core, and over the years a great deal 
has been learned about the properties of the gas compris- 
ing the broad emission line region (BELR) . It is photoion- 
ized, since reverberation studies (e.g., Clavel et al. ll991|l 
show the direct response of emission line strengths to 
continuum variability. The absence of deep Lya absorp- 
tion indicates that the BELR covers only 5 — 25% of 
the continuum source (e.g., Bottorff et al. ll997|l . It has a 
small volume filling factor ~ 10~^ (as determined from 
the observed line strength to continuum ratio; Netzer 
I199U|I . It also generates a wide range of line profile shapes 
(indicating that the geometry and kinematics are com- 
plex and varied), and shows evidence of at least a two- 
component structure (CoUin-Souffrin et al. ll982[ 119861 
Wills et al. ll985(l . One of these components is associated 
with high ionization lines, including C ill], C IV, and other 
multiply ionized species, and is known as the HIL. The 
second component can be identified with the low ioniza- 
tion lines which include the bulk of the Balmer lines, and 
lines of singly ionized species (e.g., Mgii, Cii, Fell) and 
is known as the LIL. 

The regions emitting the LIL and HIL display differ- 
ent kinematics, as deduced from studies of the profiles 
and Hne widths (e.g., Gaskell UnHHI Sulentic et al. lTMSji . 
and the HIL are systematically blue-shifted with respect 
to the LIL (see, e.g., Sulentic et al.UOOO). To account for 
the variability of the low ionization Mgn and Balmer lines 
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the LILs must be optically thick (e.g., Ferland et al.! 1992ll . 
On the other hand, optically thin gas may account for the 
Baldwin effect (a negative correlation between the ultra- 
violet emission-line equivalent width and continuum lumi- 
nosity), although it remains to be seen if this is due to sam- 
ple biases (Sulentic et al. l2U(JU|l . and for the Wamsteker- 
Colina effect (a negative correlation between between the 
Civ A1549/Lya ratio and continuum luminosity; Shields 

et ainnnsii- 

It is now fairly clear that the Balmer lines form at 
the surface of an accretion disk, or close to the surface 
in an accretion disk wind (e.g., Collin- Souffrin et al. ll988l 
Marziani et al. ll996l Nicastro I2000|l . and approximately 
three quarters of the total luminosity of the broad-line 
emission is estimated to arise in the LILs (Collin-Souffrin 
et al.'T9S8") . The geometrical distribution and kinematics 
of the HIL gas is, however, much less clear: it may arise 
in a spherical outflow, but a biconical 'jet-like' distribu- 
tion is another possibility. BELR size determinations are 
commonly obtained using reverberation and photoioniza- 
tion techniques, though it was concluded in a recent study 
that gravitational microlensing could be a useful alterna- 
tive method, particularly for the HILs (Abajas et al. 2002 ). 

Many theoretical explanations have been proposed for 
the origin of the BELR. They include: i) magnetic acceler- 
ation of clouds off accretion discs (Emmering et a\. J332} \ 
ii) the interaction of an outflowing wind with the surface 
of an accretion disc (Cassidy & Raine '19961 ; iii) interac- 
tion of stars with accretion discs (Zurek et al. 1994j; iv) 
tidal disruption of stars in the gravitational field of the BH 
fRoos ni^92|l : v) interaction of an AGN wind with super- 
novae and star clusters (Perry & Dyson 1985; WiUiams & 
Perrv [T^M|l : vi) emission from accretion shocks (Fromerth 
& Melia l^UCT l : vii) ionized stellar envelopes (e.g.Torricelli- 
Ciamponi & Pietrini 2002J . Many other models have been 
shown to possess serious difficulties (see references in 
Pittard et al. l2()()l|l : in particular, any model must over- 
come the 'confinement problem', and/or continually gen- 
erate clouds. 

Many of the mechanisms on which the various models 
are based will probably contribute to the production of the 
observed BELR gas. However, it is clear that some of the 
proposed mechanisms will be more dominant than others, 
at least under certain conditions. For example, the rate of 
tidally disrupted stars in high luminosity AGNs is likely 
to be too low to account for much of the BELR in these 
objects. Which are generally the dominant contributions 
remains, to date, largely unknown. 

Even if the SNR- AGN wind interaction is not the dom- 
inant formation mechanism of the HIL gas, it is of inter- 
est due to its potential as a diagnostic of the AGN wind. 
Supernovae must occur close to the central AGN engine 
- in our own Galaxy there exists a cluster of a few dozen 
evolved massive stars with initial masses M > 20 M0 (see 
Figer & Kim 120021 and references therein) in a region of 
1.6 pc diameter centered on Sgr A*. Recent absorption line 
measurements from one of the high velocity stars are con- 
sistent with an 08-BG dwarf with a mass ~ ISM© and 



a highly eccentric orbit which brings it within 1900 AU 
(w 6.3 X lO^'' pc) of the supermassive black hole (Ghez et 
al. l2003|l . We can expect a similar, if not more extreme, 
situation in the central regions of AGN. 

The evolution of SNRs in a high density static ambient 
medium has been studied by Terlevich et al. I|1992|l . with 
particular application to the formation of BELRs in star- 
burst models developed to obviate the existence of super- 
massive black holes in the centres of AGNs. More recently, 
the additional influence of an intense continuum radiation 
field on the evolution of SNRs has been examined (Pittard 
et al. 12001(1 . With Compton cooling and heating processes 
included in these calculations, the powerful flux of ioniz- 
ing radiation influences the thermal evolution of shocked 
regions. A central flnding was that shocked gas could radi- 
ate efficiently enough to cool to temperatures and densities 
appropriate for the HIL. 

In this paper we present calculations which extend the 
work of Pittard et al. (7UUT). We describe the results of 
2D axisymmetric hydrodynamical models of the interac- 
tion of an AGN wind with a supernova remnant. As the 
formation, evolution, and structure of cold gas is of par- 
ticular interest, we have determined the mass of cool gas 
as a function of time and present some simple modelling 
of line proflles. In a future paper, we will perform more 
detailed line proflle modelling and will compare the results 
closely to observations. 

In Sec. 12 we discuss the details of our calculations; in 
Sec. 13 we discuss our results; in Sec. 0] we present some 
preliminary line profile modelling; and in Sec. El we sum- 
marize and discuss future work. 

2. Details of the calculations 

We performed our calculations with an adaptive grid hy- 
drodynamical code which is second order accurate in space 
and time (see, e.g., Falle & Komissarov 119961 11998|l . We 
used a Chevalier-Nadyozhin similarity solution specified 
by n = 12, (5 = 0^ as the initial profile for the SNR 
f Chevalier 119821 Nadvozhin 1 1 985|l . In the majority of our 

^5 = specifies expansion into a constant density ambient 
medium. In reality the ejecta initially expands into circum- 
stellar gas expelled from the progenitor by its wind and by 
any outbursts. The high ambient pressures considered in this 
work will confine this material to a smaller volume than is the 
case for massive stars in a typical Galactic environment. For 
instance, with canonical values of M = lO~^M0yr~^ , and 
V — 2000 km s~^ for the progenitor wind, and n = 10^ cm~^ 
and T = 10^ K for the surrounding environment, the ram pres- 
sure of the progenitor wind is balanced by the thermal pressure 
of the surrounding environment at a distance of 8.5 x 10^* cm 
(2.75 X 10^* pc) from the center of the progenitor. If the am- 
bient medium flows past the progenitor, substantially smaller 
coneflnement radii can occur on the upstream side. For the 
same parameters as above and v = 3000 km s~^ for the AGN 
wind, ram pressure balance occurs at a distance of 10^* cm 
(3.2 X 10~^ pc) from the center of the progenitor. These dis- 
tances are much smaller than the size of the subsequent rem- 
nants (see Sec.|HJ, thereby justifying our assumption of a con- 
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calculations we assumed a canonical explosion energy of 
10^^ ergs and ejecta mass of 10 M©, which is typical of a 
type II SN. 75 per cent of the mass and 58 per cent of 
the explosion energy are contained within a constant den- 
sity core. Heating and cooling rates for a canonical AGN 
spectrum were kindly supplied by Tod Woods (cf. Woods 
et al. ll996|l and are included in our calculations. They are 
valid in the optically thin, low-density regime with solar 
abundances. Further details and assumptions can be found 
in Pittard et al. ^2001") and references therein. 

The thermal equilibrium of gas irradiated by the in- 
tense continuum in an AGN may be described in terms of 
several ionization parameters. When cold and hot phases 
exist with comparable pressure, it is convenient to use a 
definition based on the ratio of the ionizing photon pres- 
sure to gas pressure (cf. Krolik, McKee & Tarter [BjS 1(1 : 

where r is the distance to the central continuum source, 
n is the gas density, T the gas temperature, Li, the dif- 
ferential luminosity, j^l the frequency at the Lyman limit, 
and the other symbols have their usual meanings. Another 
common definition is the ratio of the ionizing photon den- 
sity to the gas density. 

The relation between these two ionization parameters 
depends on the shape of the ionizing spectrum. For 
the canonical AGN spectrum in Cloudy^ (see Woods et 
al. ll996) . which we adopt for this work, 

U = 0.0182r4S, (3) 

where T4 is the gas temperature in units of lO^K. 
Comparison with Roos H1992|l shows that our adopted 
spectrum is neither particularly soft or hard. 

In Figure ^ we show the thermal equilibrium curve for 
the assumed AGN spectrum as a function of temperature 
and ionization parameters, S and U. At low temperatures 
photoionization heating and cooling due to line excitation 
and recombination are in near balance. At high tempera- 
tures, the equilibrium arises from the balance of Compton 
heating and cooling. The exact shape of the thermal equi- 
librium curve at intermediate temperatures is a compli- 
cated function of the irradiating spectrum, the assumed 
abundances and thermal processes (cf. Krolik et al. .l981J . 
and varies substantially from source to source. Since we 
are not modelling a specific object, we do not concern 
ourselves with the details of this part of the equilibrium 
curve. 

stant density surrounding medium and 5 = 0. We note, how- 
ever, that a more thorough treatment would consider the eflFect 
of the mass of the confined wind on the results presented in 
Sec. 13 since this could be significant compared to the amount 
of mass in the surrounding medium which is swept up by the 
time that the ejecta core interacts with the reverse shock. 
^ Per land EUnri 
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Fig. 1. Thermal equilibrium curve for the standard AGN 
spectrum in Cloudy (see Woods et al. ll99fi|l . 



To obtain cool gas in thermal equilibrium we require 
ionization parameters S < 10. As noted by Perry & Dyson 
( 1985) , shocked gas cooled back to equilibrium can have 
a value of S much lower than its pre-shock value. This is 
because the post-shock density and pressure can be much 
greater than the pre-shock value. Therefore, strong shocks 
can create conditions for the gas to cool to temperatures 
much lower than the surrounding ambient temperature. 
The crucial question is whether the shocked gas remains 
at high densities and pressures for long enough to cool 
from its post-shock temperature to T ~ 2 x 10^ K. In our 
earlier ID work (Pittard et al. i2001|l we demonstrated that 
this was possible for a supernova in a characteristic AGN 
environment. 

In the next section we present results from 2D 
simulations of a SNR evolving in an AGN environ- 
ment. We have computed models with ambient densi- 
ties n = lO^jlO^cm"'^ and with AGN wind speeds v = 
0, 3000, 5000, 7000km s~^ . The initial radius, expansion 
speed, and age of the SNR in our models is specified in 
Table^for each of the ambient densities. All models have 
the same ionization parameter and temperature for the 
ambient medium (S ss 150, T = 1.33 x 10^ K) unless oth- 
erwise stated. We further assumed that the central con- 
tinuum source is distant enough that the flux of ionizing 
radiation is constant over our computational volume. The 
SNR was evolved until the pressure of the shocked gas 
drops to the point where it is no longer able to exist in 
the cool phase (T ~ 10^ K). As the remnant expands we 
periodically regrid our model to a coarser set of grids. 
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Table 1. Initial radius, R, expansion speed, Wcxp, and age, 
t, of the SNR for models with E = 10^^ erg and M = 
10 Mq, as a function of the ambient density, n. 



n(cm 


RilQ--' pc) 


t;cxp(kms~^ ) 


t (yr) 


10'' 


3.4 


11,000 


0.3 


10^ 


1.2 


12,000 


0.1 



3. Results 

3.1. Expansion into a stationary medium 

In Figs. 121 and 131 we show the evolution of a SNR expand- 
ing into a stationary environment with n ~ 10^ cm~'^ . 
The shocked gas rapidly loses energy, first through inverse 
Compton scattering, and then through free-free and line 
cooling, and is compressed into a relatively thin zone with 
the unshocked ejecta dominating the remnant volume. 
The shocked gas has a fractional thickness of 0.132 in the 
adiabatic self-similar solution ('Chevalier ll982l Nadyozhin 
11985(1 . but at t = 5 yr the radiative energy loss has reduced 
this to just 0.036. 

Since the shocked ejecta are denser than the swept-up 
gas, they cool quicker and form a thin dense shell within 
the narrow region of shocked gas, which is bounded on 
its interior surface by the reverse shock. Since the dense 
shell is decelerated by the hot gas on its outside surface, 
it is subject to the Rayleigh- Taylor instability, and there 
have been many investigations of this behaviour in SNRs 
(e.g., Chevalier et al.[TW, Chevaher & BlondinTMB; Jun 
et al. lTMSI Blondin et al. 2001, Blondin & Ellison 2001 
Wang & Chevalier l^nna . 

While previous numerical calculations of remnants in 
the literature showed that the shocked ejecta were un- 
able to distort the position of the forward shock (the 
Rayleigh- Taylor (R-T) "fingers" being limited to about 
half the thickness of the region of hot, swept-up gas; 
Chevalier et al. Ijl992fl . a number of exceptions have re- 
cently been found. For instance, the existence of circum- 
stellar cloudlets was found to enhance the growth of the 
R-T fingers by generating vortices in the swept-up gas 
(Jun et al. ll996|l . Alternatively, if the forward shock is an 
efficient site for particle acceleration, the shock compres- 
sion ratio is increased and the region of hot, swept-up gas 
is reduced in thickness, allowing the convective instabili- 
ties to reach all the way to the forward shock (Blondin & 
Ellison l?UU 1(1 . Finally, if the ejecta are clumpy, the greater 
momentum of the clumps enables them to push through 
the position of equlibrium pressure balance and to perturb 
(or puncture) the forward shock (cf. Blondin et al. l2001l 
Wang & Chevalier EOnSl- 

To this list we can add the present work. In the high 
density, high radiation flux environment which we con- 
sider, the region of swept-up gas is significantly more ef- 
ficient at radiating energy than in less extreme environ- 
ments, and is thus more strongly compressed than nor- 
mal. In this sense our models mimic the higher compres- 



sion ratios found in models with efficient particle accelera- 
tion (Blondin & Ellison .2001.) , and like them we find that 
the high density "fingers" are able to distort the forward 
shock. Chevalier et al. H1992fl had previously suggested the 
possibility of a highly radiative inner shock front to ex- 
plain protrusions seen in VLBI observations of SN 1986J. 
The large decrease in entropy of the shocked ejecta greatly 
enhances this instability. 

At < = 10 yr, the R-T "fingers" have grown so long, 
and distorted the forward shock to such an extent, that 
the instability begins to resemble the non-linear thin shell 
instability (hereafter NTSI; Vishniac ri994fl . To our knowl- 
edge this has never been seen before in SNRs, although 
it is a common phenomena in simulations of wind blown 
bubbles (e.g., Garcia-Segura et al. ll996|l and coUiding stel- 
lar winds (e.g., Stevens et al. ll992l Pittard et ai ling . 

We note that the "fingers" do not form at constant in- 
tervals along the thin shell. Since they are not deliberately 
seeded, they naturally develop from noise within the code 
itself. It is known from both analytical and numerical work 
that small wavelength modes grow most rapidly for both 
the R-T instability (Chandrasekhar ITMTl Youngs ITUMjl 
and the NTSI fVishniac [TUMI Blondin & Marks [TMB|l . 
However, it is extremely diflacult to resolve modes of the 
thin-shell instability, as a sufficient number of grid cells 
across the thin shell is needed to follow the tangential 
flow of material (Mac Low & Norman 1993 ). Given that 
this is not the case in our models, we expect the develop- 
ment of this instability to differ in higher resolution runs. 
Nevertheless, we do not expect it to drastically alter the 
mass of gas that cools (see Section j^J. 

Until about t = 7yr, the ejecta envelope impacts the 
reverse shock, and the pre-shock density remains relatively 
constant at n « 7 x 10^ cm^'^ . However, soon after this the 
reverse shock interacts with the ejecta core. During this 
stage the ram pressure of the ejecta on the shocked gas 
rapidly declines as the pre-shock density now decreases as 
t~^, and the shocked region depressurizes as the reverse 
shock travels towards the center of the remnant. This loss 
in pressure raises the equilibrium ionization parameter, S, 
of the shocked gas to the point that the equilibrium tem- 
perature changes from ^ 10^ K to ^ 10^ K. Gas that has 
managed to cool to T 10"* K is then heated. This be- 
haviour is seen in Fig. O At t = 5 yr much of the shocked 
ejecta has cooled to T « 2 x 10* K, and at t = 10 yr, 
this has significant structure from the action of the in- 
stabilities. However, hy t = 20 yr most of the cool gas 
has disappeared, with only a few regions concentrated at 
the densest points of the remnant, these typically being 
at the extremeties resulting from the "fingers" mentioned 
earlier. By i = 40yr even these regions have disappeared, 
and the remnant is at an almost uniform temperature 
T « 1.3 X 10^ K. 

Remnants expanding into lower density surroundings 
are characterized by a thicker region of swept up material; 
those expanding into higher density have a thinner region 
of swept up material. Fig. ^ shows the morphology of a 
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remnant expanding into a stationary medium with n = 
lO^cm-^ . 

3.2. Expansion into an AGN wind 

Several major differences in the morphology and evolution 
of the remnant are seen when it expands into a high ve- 
locity flow rather than into a stationary medium. Fig. \^ 
shows density and temperature plots of a remnant expand- 
ing into a flow with n — 10^ cm^^ and v = 3000 km s^^ . 
Since the expansion velocity of the remnant is intially 
much higher than this, the flow velocity of the surround- 
ing medium has little effect on the remnant morphology 
at early times. However, as the expansion velocity of the 
remnant slows, the flow velocity becomes increasingly sig- 
nificant, and the remnant is distorted from a spherical 
shape in the way that one would expect. The resulting 
morphology can be compared to that seen in adiabatic 
remnants expanding into a plane-stratified (stationary) 
medium (Arthur & Falle ll995|l . The models show increas- 
ing distortion of the remnant with wind speed. 

Motion of the ambient medium also affects the devel- 
opment of hydrodynamic instabilities. On the leading edge 
of the remnant, instabilities are more vigorous, since the 
shocked region is compressed to higher densities. On the 
trailing edge, the shocks are much weaker, and the shocked 
gas is both less dense and broader in extent, which severely 
suppresses the activity of the instabilities in this region. 
Since the trailing edge has relatively low pressure, S is 
much higher here than at the leading edge, and cool gas 
forms preferentially in the upstream direction. We find 
that the evolution of the mass of cool gas is surprisingly 
similiar over a wide range of ambient densities and flow 
speeds (see Sec. 13. 4|) . 

3.3. Individual cool clumps 

In the upper panels of Fig. El we zoom in on the remnant 
limb from the simulation at t = 10 yr as shown in Figs. El 
and O The forward shock bounds the right side of the 
shocked gas, and the reverse shock bounds the left. The 
high densities apparent in the upper panels lead to rapid 
cooling of the post-shock flow, such that the morphology 
resembles that resulting from the ram-ram pressure in- 
stability. The lower panels of Fig. j^l display the remnant 
limb at t — 30 yr for expansion into an ambient density 
n ~ 10^ cm~^ (see Fig. 0}. Compared to the situation 
shown in the upper panels the post-shock flow cools far 
less strongly, and we see that the cool regions are embed- 
ded in a much thicker region of shocked gas. 

Of note is the fact that the cool regions (which we shall 
identify henceforth as "clouds" ) contain within them a 
range of densities, equilibrium ionization parameters, and 
velocities. The clouds, being surrounded by a confining 
medium, are indeed long-lived: they remain as distinct 
entities until their pressure drops to the point where their 
equilibrium temperature corresponds to the hot phase. 



Within an AGN as a whole where many young SNRs will 
exist at a given time, cool clouds will be continuously cre- 
ated (through gas cooling in the younger remnants) and 
destroyed (through gas re- heating in the older remnants). 

3.4. Evolution of cool gas 

The mass of the cool clouds as a function of time is 
shown in Figs. □ and |H1 and Table [3 We define M = 
/|^QM(S)dS. The first gas to cool below T = 10*^ K is 
initially out of pressure equilibrium with its surround- 
ings as the cooling timescale of this gas is shorter than 
the dynamical timescale of the surrounding material. This 
can be seen by examining Fig. |3 where during the re- 
establishment of pressure equilibrium, S decreases. For the 
remaining evolution, continued expansion of the remnant 
causes S to increase towards the ambient value (S « 150). 
While S is below the value separating the cool and hot 
phases (S « 30 - see Fig. ^ , gas continues to cool and the 
total mass of cool gas continues to rise. However, when 
S > 30, the clouds are subject to net heating, and the mass 
of cold gas decreases until the clouds are eventually com- 
pletely destroyed. Fig. (71 shows that not all of the clouds 
transit from a cool to a hot phase at exactly the same 
time. Pressure and density enhancements in the remnant 
caused by the action of instabilities are able to maintain 
some clouds in their cool phase for a longer duration. 

When the surrounding medium is in motion relative to 
the progenitor of the explosion, the increased compression 
of the leading edge of the remnant results in cool gas form- 
ing more rapidly and with a lower value of S. The cool gas 
is also able to survive for a longer time in such cases. For 
n — 10^ cm~^ , the total mass of cool gas is fairly insensi- 
tive to the flow speed of the ambient medium, peaking at 
M ^ 2 — 4Mq approximately 15 — 20 yr after the explo- 
sion. On the other hand, the increased compression caused 
by a wind can significantly enhance the mass of cool gas 
in situations where its formation is otherwise marginal. 
For example, when n = 10^ cm~'^ , gas barely manages to 
cool below T « 10^ K when the medium is stationary, but 
in the model with n = lO^cm"^ and v = 3000 km s~^ 
a substantially greater mass of gas cools and exists for a 
significantly longer period. 

The results of most observational work are reported 
in terms of the ionization parameter U (Eq. EJ. Since 
we know the relationship between U and S for the AGN 
spectrum adopted in our models (Eq. jSJ, we can there- 
fore easily compare our results to observations, where a 
large range in U is seen. Low ionization lines (such as 
Mgii A2800) typically have U ^ 0.01, which translates 
into S ~ 0.3 for gas at T « 20000 K. In contrast, high 
ionization lines (such as C iv A1549) are characterized by 
larger values of U - low values of Cm A977/Ovi A1034 
imply the presence of gas with U ^ 1 (Laor et al. ll994fl . 
This translates into S ^ 25 for gas existing in the cool 
phase, which is clearly in good agreement with the results 
presented in Fig. 
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Table 2. Time evolution of the mass (Mq) of cool clouds in remnants expanding into a surrounding medium with 
various density, flow speed, and ionization parameter. Columns 2-5 show results for an ambient density n — 10^ cm^"^ . 
Columns 6-11 are for an ambient density n — 10^ cm^'^. In columns 10 and 11 we give results for stationary sur- 
roundings but with either an enhanced AGN flux (causing a lOx increase in the ionization parameter of the ambient 
medium to S « 1500), or remnant parameters suitable for a type la SN explosion (see Sees. 13.51 and 13.61 respectively! . 
Dashes indicate that no cool gas exists. The bottom row contains values for the integrated area under the mass vs. 
time curve, in units of Mq yr. 
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3.5. Effect of an increased AGN flux 

The simulations shown so far were computed with an AGN 
flux which was high enough to maintain the ambient gas 
in the hot phase (i.e. S w 150 - see the equilibrium curve 
in Fig. This resulted in a relatively low value of S 
for the shocked gas and provided a good opportunity for 
the shocked gas to cool to T « 2 x 10'' K. However, if 
a remnant were to be immersed in ambient gas with a 
higher ionization parameter, the ability for shocked gas 
to cool would be compromised. This is shown in Fig. |3 
where a remnant is expanding into a stationary medium 
with S w 1500 (i.e. S is lOx higher than in Figs. [5] and O. 
The increased photon flux from the AGN increases the 
rate of Compton heating, and leads to a reduction in the 
net cooling rate. While the radius of the remnant is largely 
unchanged, its morphology is affected: the shocked gas is 
not compressed as much, and cool clouds do not form. This 
behaviour is consistent with the earlier work presented in 
Pittard et al. plHT]! . 

3.6. Effect of the stellar envelope 

To investigate whether cool clouds could form in a rem- 
nant from a type la SN explosion, we have computed 
an additional model with appropriate parameters: E = 
10^1 erg, M = 1.4 Mq, and n = 7 (cf. Chevalier |T982) . 
The main difference with this model is that the increased 
value of E /M results in higher expansion velocities of the 
ejecta and hotter post-shock gas, and a lower equilibrium 
ionization parameter. However, the remnant also deceler- 



ates more rapidly, causing S to increase at a faster rate 
than for the n = 12 case. This ultimately hinders the for- 
mation of cool gas relative to the n = 12 case, and only a 
small amount is able to form (see Table EJ- 

Recently, it is has become clear that there exists a class 
of type II supernova explosions which are under energetic 
(e.g., Zampieri et al.)2003). For two explosions examined 
in detail, Zampieri et al. ( 2003 ) found M > 14 and 
E « 0.6 — 0.9 X 10^^ erg. We do not expect remnants with 
such parameters to evolve significantly differently to our 
canonical models with M = 10 Mq and E = 10^^ erg. 

4. Line profiles 

We have calculated synthetic line profiles for emission 
from the cooled gas. The 2D axisymmetric grid is rotated 
onto a 3D cartesian grid and the emission from volume 
elements containing cool gas was integrated under the as- 
sumption that it is optically thin and the volume emission 
rate varies as n^. Since the gas is cool, thermal Doppler 
broadening is negligible. A previous investigation (Bottorff 
et al.POOO ) failed to reach any strong conclusions concern- 
ing whether microturbulence"^ was favoured by observa- 
tions, so it is not included in our model. The emission is 
blue- or redshifted according to the line of sight velocity 
of the gas. Absorption was also assumed to be negligible. 

In Fig. 1101 we show the line profiles resulting from a 
remnant expanding into an ambient medium with 3 sets 

^ Defined as any velocity field that occurs over distances that 
are small compared to a photon's mean free path. 
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Fig. 8. The evolution of the mass of cool clouds in rem- 
nants expanding into a surrounding medium with n = 
10^ cm~^ (top) and n = 10® cm~'^ (bottom). In each panel 
the mass of cool gas is shown as a function of t for flow 
speeds of u = 0, 3000, 5000, 7000 km s'^ . 



of parameters. The normalization of the line profiles has 
been set so that the peak emission is approximately 1.0, 
and we have preserved the relative scaling between the 
models (e.g., the central intensity of the n — 10® cm"'^ , 

V = Okms~^, t = 20 yr profile is approximately 4x 
greater than the central intensity of the n ~ 10^ cm^'^ , 

V = km s^^ , t = 20 yr profile). The bottom row in Fig. 1101 
shows the line profile which results if we integrate over the 
age of the remnant. In effect we sum each of the profiles in 
the rows above with an appropriate weight which reflects 
the time between each "snapshot" . 

In general, the normalization first rises and then falls 
with time as cool clouds are created and then destroyed. 
The width of the line decreases with time as the expan- 
sion speed of the remnant slows. The line is clearly flat- 
topped, which is characteristic of emission from a geomet- 
rically and optically thin spherical shell (see, e.g., Fig. 3 
in Capriotti et al. ll980l with I'min/wmax ~ 1-0). As the 
effective thickness of the "shell" in Fig. |21 is minimal, a 
tangential line of sight does not intercept an increased 
number of clouds, and "horns" are not seen in the profile. 
In contrast, the line profiles displayed in the rightmost 
column of Fig. 1101 show an increasingly rounded or trian- 



gular profile as the remnant ages, and the time-averaged 
profile displays a distintly rounded top. In this model the 
remnant is expanding into an AGN wind and the line pro- 
files (which are sensitive to the spatial distribution of cool 
gas) refiect the increasing distortion of the remnant as it 
expands. 

It is clear from Fig. |21 that the emission is from a rela- 
tively small number of clouds (particularly at later times, 
e.g., t — 20 yr). This introduces a great deal of small scale 
structure into the line profiles which we have smoothed 
out in Fig. 1101 by averaging over many different lines of 
sight. Observed line profiles from AGN are in reality very 
smooth, and it has been concluded that the number of 
emitting clouds must be ^ 10^ (Arav et al. 11997(1 . although 
this would be reduced if there was significant microtur- 
bulence. Alternatively, electron scattering could help to 
explain smooth broad line profiles, especially in the line 
wings (Emmering et al. ll992|l . In our 2D hydrodynamical 
models the clouds are in fact rings, and numerical viscos- 
ity and the finite number of grid cells limit the number 
of distinct clouds which form. Increased numerical resolu- 
tion and 3D simulations will produce many more distinct 
clouds, but present limitations mean that the line profiles 
from our models are not as smooth as seen in observations, 
and their fine structure should be ignored. 

In Fig. 111! we show the line profiles resulting from 
a remnant expanding into an AGN wind of speed v = 
3000,5000, and 7000 kms'^ for a viewing angle, 9 = 
135°. Here 9 is defined as the angle between the AGN 
wind vector and the vector from the observer to the rem- 
nant (i.e. 9 = 0° corresponds to the observer facing the 
side of the remnant expanding into the oncoming AGN 
wind). Inspection of Fig. [S] reveals that cool clouds ex- 
ist over a wide range of angles, namely from 9 ^ 0° to 
9 « 135°. For a line of sight with 9 = 135° (i.e. where 
the flow has a velocity component towards the observer), 
the blue wing of the line shows the greater extension at 
t = 10 yr. This is expected since the trailing edge of the 
remnant is not decelerated as rapidly as the leading edge. 
However, at i = 5 yr, the cold clouds exist predominantly 
on the leading edge (the density and cooling rate are high- 
est here) and it is the red edge of the line profile that is the 
most extended. This behaviour is also seen when the rem- 
nant expands into a flow with v = 5000 or 7000 km s^^ . 
Hence as the remnant expands, the resulting line profile 
flips from a red to a blue-shift (and vice- versa for 9 < 45°). 

For each of the 3 distinct cases of ambient medium 
investigated in Fig. 1101 the time-averaged line profile 
displayed in the bottom row has a FWHM w 4500 — 
5000 km s~^, and a full width at zero intensity of 
6000 — 8000 km s^^ . This range lies almost in the mid- 
dle of the distribution of FWHM measures of the broad 
component of Civ A1549 for radio-quiet sources (Sulentic 
et al. l^ITO . 

In Table O we list various statistics for the time- 
averaged line profiles shown in Fig. ^2 Radio quiet sources 
show a scatter in the asymmetry values for the broad 
component of the CivA1549, with A.I. = ±0.2, and 
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Table 3. Statistics for the time-averaged line profiles 
shown in Fig. ^Jl"- — 10^ cni~^ ). Detailed are the profile 
centroid measures at different levels of the line peak, the 
full- width at half maximum (FWHM) , and the asymme- 
try index, A.I. = (C(|) - C(i))/FWHM. The columns 
labelled wsooo, I'sooo, and vtooo refer to the models where 
the AGN wind speed is 3000,5000 and 7000 km s'^ re- 
spectively. 



Statistic 


^'3000 


«5000 


U7000 


C(l) (kms-^) 





900 


1000 


C(3/4) (kms-i ) 





1400 


1200 


C(l/4) (kms"' ) 


-100 


500 


1200 


FWHM (kms-^) 


4200 


3200 


4200 


A.I. 


0.0 


0.3 


0.0 



(A.I.) « +0.04 (Marziani et al. HlMjl . This lack of a strong 
preference for red or blue asymmetries is currently much 
in need of confirmation, especially when one considers the 
systematic blueshift that C ivA1549 shows in RQ sources 
(see Sulentic et al. l2()()()l and references therein). It is in- 
teresting to note that our models produce values for the 
asymmetry index which are compatible with the observa- 
tionally deteremined values. 

Observations have revealed great diversity in the pro- 
files of particular broad emission lines, such as C iv A1549, 
whereas variations in other lines, such as Oiv] A 1402, are 
much reduced in comparison (see, e.g., Francis et al. 1992 
Wills et ai niMjl . The variation in C iv A1549 from object 
to object can be explained if a two-component structure 
is invoked - a broad "base" which is relatively constant 
between objects, and a narrow "core" whose contribution 
to the overall line strength is variable. With this model, 
lines with large equivalent widths tend to be narrow (or 
"cuspy") and have large peak fluxes, as observations re- 
quire. In contrast, lines such as Oiv] A1402, which show 
far less variability from object to object, are essentially 
composed of a single broad component. In such cases, the 
contribution from a variable narrow core is small. 

The observational interpretation of broad line profiles 
has advanced greatly over the last decade, partly due to 
a better understanding of sample biases and more careful 
consideration of contaminating lines and superposed nar- 
row components. Sometimes an obvious inflection demon- 
strates a superposed narrow component. Alternatively, if 
a line has a rounded top, one can infer that the narrow 
component is largely absent. 

Inflections also support the interpretation that the 
broad lines are formed in several kinematically and/or ge- 
ometrically distinct emitting regions, as does the frequent 
"mismatch" in profile wings (Romano et al. [T^9 6). One 
possibility is that the variable narrow component forms on 
the outer edge of the BELR (this is sometimes referred to 
as an "intermediate line region"), with the broad compo- 
nent formed at smaller radii within the BELR. Variations 
are then most likely explained as a changing covering fac- 



tor (if optically thick) or volume cmissivity at the outer 
edge. Support for this model comes from variability mea- 
surements of NGC 5548 and NGC 4151 CClavel lTMTbll . in 
which the line cores (defined as the central 3000 km s~^ ) 
lag continuum changes with approximately twice the delay 
of the wings. An alternative model consisting of a biconical 
outflow and an accretion disc continuum source has been 
shown to have problems in explaining the near-constant 
velocity width and equivalent width of the emission line 
wings (Francis et al. ll992|l . and there are doubts concern- 
ing the dependence of the line equivalent width on orien- 
tation via. an axisymmetric continuum, which this model 
requires (Wills et al. ITU^ . 

5. Discussion 

The results in Sec.|31demonstrate that it is possible to cool 
shocked supernova ejecta to T 10"* K in the inner regions 
of a QSO. Although our results differ from the original 
proposals of Dyson & Perry (1982) and Perry & Dyson 
JEHHIl, which were for the shocked ambient medium to 
cool, the resulting cool gas nevertheless has properties 
(densities, column densities, velocities and ionization pa- 
rameters) compatible with those inferred for gas emitting 
the high ionization lines in QSOs. 

A parameter space study shows that for ambient den- 
sities of n = 10''' — 10^ cm^^ , about 1 — 4Mo of material 
can cool to low temperatures. Such gas then persists for 
~ 10 — 20 yr before the continuing expansion of the rem- 
nant reduces the density and pressure of the cool gas to 
the point where its equilibrium temperature and ioniza- 
tion parameter corresponds to the hot phase. This result 
is robust for a wide range of velocities of the surrounding 
medium, although the spatial distribution of the cool gas 
around the limb of the remnant, and hence the resulting 
optical/UV line profile, is sensitive to this detail. We find 
that the integrated value of the mass of cool gas over the 
remnant lifetime generally increases with the density and 
the flow velocity of the surrounding medium. The highest 
value found in our simulations, Mt = 77MQyr, is ob- 
tained when n = 10^ cm~^ and v ~ 5000 km s~^ . 

A supernova rate of 1 /yr would then imply a mass for 
the clouds emitting the HILs of up to ~ 80 M0 , This is eas- 
ily compatible with the lower end of BELR mass estimates 
in the literature (e.g., Peterson ll997|l . although our model 
(and most others) would be severely challenged to explain 
much more extreme estimates of the mass of BELR gas 
(see Baldwin et al. l20()3l and references therein). We note 
that it is currently unclear how this mass is partitioned 
between the HIL and LIL gas in these higher estimates. 

In earlier work it was shown that for typical QSO pa- 
rameters the power going into supernova remnants is com- 
parable to that of the QSO wind, but is much less than 
the bolometric QSO luminosity (Perry & Dyson 11985(1 . 
However, it is more difficult to estimate whether emission 
from the SN would be visible above the QSO in a specific 
waveband. The typical J-band magnitude for a QSO at a 
redshift z ~ 1 is ~ 18 — 19, whereas the J-l-H band mag- 
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nitude for a type la SN is ^ 24 at comparable z. On this 
basis, individual SN will not be discernible, but clearly this 
conclusion depends on the luminosity of the QSO, as well 
as other variables such as the orientation of the SN with 
respect to the molecular torus, and the ambient density of 
the surroundings (e.g., if the SNR expands into a nearby 
molecular cloud then its luminosity could be significantly 
increased). Detailed numerical modelling will be required 
to determine the likelihood of this possibility. 

One of the most interesting questions concerning 
AGNs is the connection between nuclear and starburst 
activity. Mixed starburst-AGN sources may be recogniz- 
able as outhers in the EW(H /3) - EW(Fe ii A4570) plane 
fSulentic l2000| . In this sense AGNs displaying galaxy-wide 
starburst activity are in the minority. However, to provide 
enough fuel for the high luminosity QSOs, the accreting 
medium needs to be particularly dense. Possible mecha- 
nisms to augment the density include the winds and ex- 
plosions of massive stars, and galaxy collisions. As the lat- 
ter are often associated with starbursts, vigorous massive 
star formation in the inner regions of AGN is probably a 
necessity. If a significant component of the BELR results 
from cool gas in SNRs (as explored in this paper), then 
essentially all AGNs must have a nuclear starburst. 

We have assumed in our models that both the ejecta, 
and the surrounding interstellar medium, are homoge- 
neous and have solar abundances. While the presence of 
large-scale macroscopic mixing of ejecta in core collapse 
SNe has been well established on both observational and 
theoretical grounds (see Blondin et al. l2()()ll and references 
therein), such mixing is not complete on a microscopic 
level. For example, X-ray observations of the Vela SNR 
(Aschenbach et al. ll995fl show several fragments outside 
of the general boundary, and ASCA (Tsunemi et al. ll999|l 
and XMM-Newton (Aschenbach & Mivata .2(7n5|l obser- 
vations of fragment A have revealed a significant over- 
abundance of Si and Mg, confirming that this fragment 
is ejecta. Widespread evidence that the ejecta of core col- 
lapse supernovae are clumpy is further noted in Wang & 
Chevalier H2002(l . The possibility that supernovae are ex- 
plosions of "shrapnel" which give rise to a complex outer 
boundary has been discussed by Kundt l|1988|l . Density 
enhancements in the surrounding medium may also occur. 
The interaction of ejecta with a clumpy wind has been pro- 
posed as the origin of the broad and intermediate-width 
lines in the spectrum of the peculiar SN 1988Z (Chugai & 
Danziger 1994). 

Knots of X-ray emission seen in the Tycho SNR 
(a type la explosion) also indicate clumpy ejecta, and 
abundance variations between the knots indicate that 
the mixing of the deep Fe layer changes from place to 
place (Decourchelle et al. l2()()l|l . However, on a larger 
scale a general stratification is seen, with the lighter ele- 
ments found predominantly at large radii, and vice-versa. 
Emission from Fe has the smallest radius, confirming the 
onion shell structure expected from deflagration models. 

While this work has used solar abundances and homo- 
geneous media for simplicity, it is clear that the majority 



of the cool mass will be nuclear processed material and 
that there will be localized density enhancements in the 
ejecta. Such material will cool more efficiently, and higher 
masses for the cold phase will be obtained. In this sense, 
the values in Table |21 should be viewed as lower limits. 
Future models will eventually need to treat in a realis- 
tic fashion the inhomogeneities in density and abundance 
that we know exist. 

Our investigation of the influence of an AGN environ- 
ment on the dynamics and evolution of a supernova rem- 
nant is ongoing. In the next paper in this series we will 
study the dynamical influence of the QSO radiation field. 
We also note that the combined wind from a group of 
early-type stars may provide the necessary conditions for 
the formation of cool regions. We anticipate that this sce- 
nario will be more relevant in the nuclei of Seyfert galaxies, 
since supernova explosions will evacuate all but the most 
tightly bound gas in them (Perry & Dvson [r985fl . Finally, 
it is clear from our models that while the supernova-QSO 
wind interaction is conceptually simple, the BELR is likely 
to be a very complicated region in practice. 
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Fig. 2. The evolution of a SNR expanding into a stationary environment of density n = 10^ cm~^ . Shown are logarith- 
mic plots of the density with a scale from 2.4 x 10"^ cm~^ (white) to 3.1 x 10^ cm~^ (black). The age of the remnant is 
shown in the upper left corner of each plot. Note the change in the axis scaling between the plots. The unit of length 
in these and subsequent plots is 1.1 x 10^^ cm (0.0036 pc), and the ionization parameter of the surrounding medium, 
S « 150, unless explicitly stated otherwise. 
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Fig. 4. The morphology of a SNR expanding into stationary surroundings with n = 10^ cm ^ at t = 30 yr. Left: 
logarithmic density (grayscale from 4.7 x 10''cm~^ (white) to 1.8 x lO^cm"'^ (black)). Right: logarithmic temperature 
(grayscale from 10^ K (white) to 1.58 x 10^ K (black)). The unit of length is 2.24 x 10^^ cm (0.0073 pc). 
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Fig. 5. The evolution of a SNR expanding into an AGN wind with n = 10^ cm~^ and v = 3000 km . In this figure 

the AGN wind propagates down the page and the explosion site is at position (0,0). The top row shows logarithmic 
plots of the density with a scale from 2.4 x lO^cm"^ (white) to 3.1 x 10^ cm""^ (black), while the bottom row shows 
logarithmic plots of the temperature with a scale from 10^ K (white) to 1.58 x 10^ K (black). The age of the remnant 
is shown in the upper left corner of each plot. Note the change in the axis scaling between the plots. 
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Fig. 6. The density and temperature structure of individual bowshocks and cool clumps. The top panels show the 
structure existing at t = 10 yr when a remnant expands into an ambient density n = 10^ cm~^ . The bottom panels 
show the structure existing a,t t = 30 yr when a remnant expands into an ambient density n = 10^ cm""^ . The left 
panels display density gray-scales, contours, and flow vectors, while the right panels display temperature gray-scales 
and contours. 
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Fig. 7. The evolution of cool gas in remnants expanding into a surrounding medium with n = 10® cm ^ . The 
mass of cool gas with an ionization parameter below S is shown as a function of t for flow speeds of v = 
0, 3000, 5000, 7000kms-i . 
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Fig. 10. Line profiles from cool gas formed in a SNR expanding into an AGN environment with S « 150 and integrated 
over many lines of sight (i.e. it is assumed that there is a spherical distribution of remnants in an isotropic medium and 
radiation field). The left and middle columns are for stationary environments with n = lO^cm"^ and n = lO^cm"^ , 
respectively. The right column is for an environment with n = 10^ cm"*^ and v = 3000 km s~^ . In the top two rows the 
line profile at specific remnant ages is displayed, with the age of the remnant (in years) marked next to each profile. 
The solid line in the top row corresponds to t = 20 yr for n = 10^ cm~3 ^ and t = Syr for n = 10^ cm~'^ . The remaining 
profiles in the top and middle rows correspond to the subsequent entries in Table |21 with the line width decreasing 
with remnant age. In the bottom row the average line profile integrated over the age of the remnant is displayed. Note 
the difference in the scaling of the intensity axis. 
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Fig. 11. Line profiles from cool gas formed in a SNR expanding into an AGN wind with n = 10^ cm~'^ , S 150, and 
V = 3000, 5000, or 7000 km s~^ . The line-of-sight is specified by ^ = 135° (i.e. there is a component of the AGN wind 
moving towards the observer) and the profiles are from remnants of age f = 5yr (top row), t = lOyv, t = 15 yr, and 
t = 20yr. In the bottom row the age integrated profiles are shown. 



